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Abstract. Recently A. Schrijver derived new upper bounds for binary codes 
using semidefinite programming. In this paper we adapt this approach to codes 
on the unit sphere and we compute new upper bounds for the kissing number in 
several dimensions. In particular our computations give the (known) values for 
the cases n = 3, 4, 8, 24. 



1. Introduction 

In geometry, the kissing number problem asks for the maximum number Tn 
of unit spheres that can simultaneously touch the unit sphere in n-dimensional 
Euclidean space without pairwise overlapping. The value of t„ is only known for 
n = 1, 2, 3, 4, 8, 24. While its determination for n = 1, 2 is trivial, it is not the case 
for other values of n. 

The case n = 3 was the object of a famous discussion between Isaac Newton 
and David Gregory in 1694. For a historical perspective of this discussion we 
refer to [6|. The first valid proof of the fact "rs = 12", like in the icosahedron 
configuration, was only given in 1953 by K. Schiitte and B.L. van der Waerden in 

m. 

In the seventies, P. Delsarte developed a method, initially aimed at bounding 
codes on finite fields (see HI), that yields an upper bound for r„ as a solution of 
a linear program and more generally yields an upper bound for the size of spher- 
ical codes of given minimal distance. We shall refer to this method as the LP 
method. With this method, A.M. Odlyzko and N.J.A. Sloane (IMl), and indepen- 
dently V.I. Levenshtein ( Il4ll ). proved rg = 240 and T24, = 196560 which are 
respectively the number of shortest vectors in the root lattice and in the Leech 
lattice. For other values of n, the LP method gives in many cases the best known 
upper bounds. However, for n = 3 and n = 4 it only gives the upper bounds 
T3 < 13 and T4 < 25. 
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In 2003, O.R. Musin succeeded to prove the conjectured value T4 = 24, which 
is the number of shortest vectors in the root lattice D4, with a variation of the LP 
method (see [15] and the survey fT9] of F. Pfender and G.M. Ziegler). 

To complete the picture, let us discuss uniqueness of the optimal point con- 
figurations. For dimensions 8 and 24, uniqueness was proved by E. Bannai and 
N.J.A. Sloane ([3|). Their proof exploits the fact that the LP method obtains ex- 
actly the aimed value. For dimension 3, there are infinitely many possible configu- 
rations. In the regular icosahedron configuration, the angular distances between the 
contact points are strictly greater than the required vr/S, hence these points can be 
moved around obtaining infinitely many new suitable configuration. This partially 
explains why the determination of T3 is difficult. On the contrary, uniqueness of 
the optimal configuration of points in dimension 4 is widely believed, but remains 
unproven. 

The LP method, which was established by P. Delsarte, J.M. Goethals and J.J. Sei- 
del in [9J, handles the more general problem of the determination of a bound for 
the maximal number 

A{n, e) = max{card(C) : C C 5"-^ with c • c' < cos for c,c' eC,ci^ c'} 

of points on the unit sphere with minimal angular distance 6. Such configurations 
of points, also called spherical codes with minimal angular distance 9, are of spe- 
cial interest in information theory. The kissing number problem is equivalent to the 
problem of finding yl(n, vr/S). 

In this paper, we define a semidefinite program (SDP for short) whose optimal 
solution gives an upper bound for A{n, 6) and strengthens the LP method. Com- 
putational results show that for several values of n this SDP method gives better 
upper bounds for r„ than the LP method. 

To be more precise, let us recall that the LP method relies on the existence of 
polynomials PJ^{t), satisfying the so-called positivity property: 

(1) for all finite C C ^ Pfc"(c • c') > 0. 

(c,c')ec2 

These polynomials arise as zonal spherical polynomials on the sphere, i.e. the zonal 
polynomials associated to the decomposition of the space of polynomial functions 
under the action of the orthogonal group 0(M"). 

The consideration of the action restricted to a subgroup H of 0(IR"), chosen 
to be the stabilizer group of a fixed point e G 5"*^^, leads us to some symmetric 
matrices 5^ whose coefficients are symmetric polynomials in three variables such 
that 

(2) for all finite C C ^ 5|?(c • c', c • c", c' • c") ^ 

(c,c',c")6C3 

where the sign 'V 0" stands for: "is positive semidefinite". The reason why we 
obtain matrices instead of functions comes from the fact that, in the decomposition 
of the space of polynomial functions on the sphere under the action of H, multi- 
phcities greater than 1 appear. In fact these multiplicities are exactly the sizes of 
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the corresponding matrices. From ([T]) and ^ we derive an SDP whose solution 
gives an upper bound for A{n, 9). 

Our approach adapts the method proposed by A. Schrijver in to the unit 
sphere whereas he obtains new upper bounds for binary codes from an SDP. His 
work can also be interpreted in group theoretic terms, involving the isometry group 
of the Hamming space F2 and the subgroup stabilizing (0, . . . , 0) which is the 
group of permutations of the n positions. It is very likely that many other spaces 
of interest in coding theory can be treated likewise. The case of non-binary codes 
was considered by D.C. Gijswijt, A. Schrijver and H. Tanaka in fTTjl. 

The paper is organized as follows: Section|2]reviews on the LP method. Section 
|3]introduces and calculates the semidefinite zonal matrices associated to the action 
of H and leading to the matrices S^. Section |4] defines the semidefinite program 
and its dual that estabUshes the desired bound. Section |5] discusses computational 
results. 



2. Review of the LP Method on the Unit Sphere 

We introduce the following notations. The standard inner product of the Eu- 
clidean space M" is denoted hy x ■ y. The unit sphere 

:= {x G M" : X • X = 1} 

is homogeneous under the action of the orthogonal group 0(M") = {O G M"^" : 
0*0 = /„}, where /„ denotes the identity matrix. It is moreover two-point homo- 
geneous, meaning that the orbits of 0(R") on pairs of points are characterized by 
the value of their inner product. The space of real polynomial functions of degree 
at most d on S*""^ is denoted by PoI<rf(S'"~^). It is endowed with the induced 
action of 0(M"), and equipped with the standard 0(M") -invariant inner product 

(/,£/) = — I f{x)g{x)duJn{x), 

where tOn = Y{n/2) "-^^ surface area of 5""^^ for the standard measure diOn- It is 
a well-known fact (see e.g. iMl Ch. 9.2]) that under the action of 0(R") 

(3) Pol<d(5"-i) = HS LH^L...!. H2, 

where H'^ is isomorphic to the 0(M") -irreducible space 

Harm^ = |/ £ "^[xi, . . . , x„] : / homogeneous, deg f = k,^ -^f = o| 

i=l * 

of harmonic polynomials in n variables which are homogeneous and have degree 
k. We set hi := dim(Harm^) = - nif)- 

A certain family of orthogonal polynomials is associated to the unit sphere. They 
will be denoted by P^, with the convention that PJ^ has degree k and is normalized 
by PJ^{1) = 1. For n > 3 these polynomials are up to multiplicative constants 
Gegenbauer polynomials with parameter A = n/2 — 1. So they are given 
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by P;^{t) = C^/^~^(t)/C^/^~\l), and the Gegenbauer polynomials can be 
inductively defined by C^{t) = 1, C^{t) = 2Xt, and 

kC^{t) = 2{k + X- l)tC^_^{t) -{k + 2X- 2)C^_2(t), for k>2. 

They are orthogonal with respect to the weight function (1 — t2-)A-i/2 j-j^g inter- 
val [—1,1]. For n = 2 the polynomials coincide with the Chebyshev polynomi- 
als of the first kind which can be inductively defined by To(t) = 1, Ti{t) = t, 
and 

Tk{t) = -2tTk-i{t) + Tk-2{t), for k>2, 

and they are orthogonal with respect to the weight function (1 — on the 

interval [—1, 1]. 

The polynomials PJ^{t) are related to the decomposition ^ by the so-called 
addition formula (see e.g. O] Ch. 9.6]): for any orthonormal basis (ei, . . . , e/i^) of 
HJ^ and for any pair of points x, y G 5*^"^ we have 



(4) Pk{x-y) = -r^^e,{x)eiiy). 

From the addition formula (01), the positivity property ([T]) becomes obvious: 

h" 

E e^(c)e^(c') = /^E E^^(^) 

k i=l (c,c')eC2 ^ i=l \ceC / 



hi 

{c,c')£C^ (c,c')GC2 k j=i 



{c,c')ec2 

Now we introduce the unknowns of the LP to be considered. For a spherical 
code C we define the two point distance distribution 

'■= 77777 card{(c, c') G : c • c' = u}, 

card(G) 

where u G [—1, 1]. Clearly, only a finite number of 2;(ti)'s are not equal to zero, 
and the positivity property can be rewritten as a linear inequality in the a;(u)'s: 

(5) E ^M^rM>0. 

«e[-i,i] 

Moreover, the number of elements of C is given by card(C) = J2ue[-i i] ^(^)- 
Noticing the obvious conditions x(l) = 1, x{u) > 0, and x{u) = for cos 9 < 
ti < 1 if the minimal angular distance of C is 9, we are led to consider the following 
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linear program: For any d > 1, the optimal solution of the linear program 

max { 1 + x(u) : 

l.cos 8] 

x{u) = for all but finitely many u G [—1, cos 6], 
x{u) > for all u G [—1, cos 9], 

1 + E„e[-i,cose] ^iu)Pki^) > for all A: = 1, . . . , d}, 
gives an upper bound for A{n, 9). The dual linear problem is 

d 



mill { 1 + ^ A 



(7) fc=i 

/fc > for all = 1, . . . , d, 

ELi fkPkiu) < -1 for all n G [-1, cos^]}. 

By the duality theorem (of. fW\) any feasible solution of (|7]) gives an upper bound 
for the optimal solution of The dual linear program can be restated in the 
following way involving polynomials: 

Theorem 2.1. (see e.g. |9, Th. 4.3], 03, US, El Ch. 9]j 

Let F{t) = Efe=o fkPk{t) be a polynomial of degree at most d in M[t]. If 

(a) fk > Ofor all k > 1 and /o > and 

(b) F{u) < Ofor all ue [-l,cos6'], 
then 

F(l) 
Jo 



3. SEMIDEFINITE ZONAL MATRICES 

Now we fix a point e G and let H := Stab(0(M"), e) be the stabilizer of 

e in 0(M"). Obviously, H ~ 0(M""i) since 0(M""i) can be identified with the 
orthogonal group of the orthogonal complement of Me. 

It is a classical result (see e.g. [24, Ch. 9.2]) that for the restricted action to H 
the decomposition of Harm^ into ff-irreducible subspaces is given by: 

k 

Harm^ ~ QHarm^'-^ 

i=0 

Hence, each of the in ^ decomposes likewise: 

where i/f^ ^ ~ Harm['"\ We give an explicit description of this decomposition 
in the proof of Theorem |3.2[ 

We summarize the situation in the following picture. 
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Pol<d(5"-i) = ± Hf ± ... ± H2 





Hf 




. ± 




Tjn—l 








Tjn—l 
^1,1 







Tjn—l I ijn—l I I un—l 

^0,0 -L -^0,1 -L ••• -L ^o,d 

Tjn-l 



^ HZ/ 

The isotypic components of the ff-decomposition of Fol<,d{S"'^^) are 

(9) Jfc:=/f^^^i±...±i7fc"/~((i-fc + l)Harm^~\ for k = 0, . . . ,d. 

Now we show how to associate to each Ik a "zonal matrix" in view of an analogue 
of the addition formula (|4]l. 

Theorem 3.1. Let 1 = Rq ± Ri ± . . . 1. Rm ~ (m + 1)R be an isotypic 
component o/Pol<rf(S'"~^) under the action of H, with R an H -irreducible space 
of dimension h. Let (eo.i, • • • , eo,/i) be an orthonormal basis of Rq and let (pi : 
Rq — > Ri be H -isomorphisms preserving the inner product on Pol<,ii{S"'~^). Let 
Cij = (pi^SQj), so that (cj^i, . . . , Ci^h) is an orthonormal basis of R^. Define 

I eo,i(x) . . . eo,/i(x) 



and 



\em,i(2;) ... em,h{x)^ 



Z{x,y) := E{x)E{yY e M{™+i)x(m+i)_ 

Then the following properties hold for the matrix Z: 

(a) Z(x, y) does not depend on the choice of the orthonormal basis of Rq. 

(b) The change of (pi to —(pi for some i or the choice of another decomposi- 
tion of I as a sum of m + 1 orthogonal H-submodules changes Z(x, y) to some 
OZ{x,y)0^ with O £ 0(M'^+^). 

(c) For ally e H, Z{g{x),g{y)) = Z{x,y). 

(d) (Matrix-type positivity property) 

(10) For all finite C C 5"-\ Z{c,c) h 0. 

(c,c')GC2 

Proof (a) If (eo,i, . . . , eo,/i) is another orthonormal basis of Rq, then there is an 
orthogonal h x h matrix O with (eo,i, • • • , eo,h) = {^o,i, • • • > eo,/i)0. In this case 
the matrix E{x) is changed to E{x)0 and, since OO* = Ih, the matrix Z{x, y) 
stays unchanged. 

(b) By Schur's Lemma and by the irreducibility of R, there are only two possible 
choices for (pi, namely (pi and —(pi, once the subspaces Ri are fixed. 

Let T = 5o -L . . . -L Sm be another decomposition of T, together with H- 
isomorphisms ipi : Ri ^ Si preserving the inner product on Pol<^(S'"^^). Then 
ip = (V'O; • • • ) V'm) defines an ff-endomorphism of I. Again by Schur's Lemma, 
for a suitable choice of basis in Ri and by permuting rows and columns, the matrix 
of Ip is block diagonal with h blocks of size (m + 1) x (m, + 1) and with the same 
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(m + 1) X (m + 1) matrix O G 0(M'"^^) as blocks. This means that E{x) changes 
to OE{x) and so Z{x, y) becomes OZ{x, y)0^. 

(c) Since eij{g~^{x)) = {geij){x), the computation of Z {g~^ (x) , (y)) 
amounts to replace in the definition of Z{x, y) the j by geij. Since Ri is H- 
stable, ejj := geij, with j = 1, . . . , /i, is another orthonormal basis of Ri, and 

= (piigeij) = g4>i{ei,j) = gcij = eij. 

Hence from (a) we conclude Z {g^^ (x) , g^^ (y)) = Z{x,y). 

(d) We have = ( Z ^(^)) ( E ^(c)) \ 0. □ 

(c,c')GC2 cGC cGC 

The orbits of H on pairs of points on the unit sphere x,y G S"^^ are char- 
acterized by the values of the three inner products e ■ x, e ■ y, and x ■ y. By 
definition the coefficients Zij{x,y) of Z{x,y) are polynomials in the variables 
xi, . . . ,Xn,yi, . . . ,yn- Then, property (c) of Theorem l3. 1 l implies that they can be 
expressed as polynomials in the three variables u = e ■ x, v = e ■ y, and t = x • y. 

By Z^, for < /c < d, let us denote the matrix associated to as defined 
above, and more precisely to the decomposition (|9l) of T^.. Now we shall calculate 
the matrix Y^{u, v, t) with 

(11) Zl{x,y) = Y^{e-x,e-y,x-y). 

Theorem 3.2. With the above notations, we have, for all < i,j < d — k, 

(12) {Y,^)^^^{u,v,t) = \i,,Pp+'\u)P;^+'\v)Ql-Hu,v,t), 
where 



^ ^^{l-v?){l-v^) 



and 



Proof. We explicitly use an orthonormal basis of H^~^]^^ to calculate yj!^(n, t). 
Such a basis is constructed in [L Ch. 9.8]. Let us recall the construction. For 
X £ let 

x = ue + \/l — n^C) 
where u = x-e and C belongs to the unit sphere 5"^^ of (Me)-^. With / G C 
Pol<fc(5'^-2) we associate (/?(/) G Pol<fe(S'"-i) defined by: 

V.(/)(x) = (l-uY/2^(C). 

Note that the multiplication by (1 — v?)^!"^ forces to be a polynomial function 
in the coordinates of x. Clearly Lp commutes with the action of H. Hence (f{Hj^~^) 
is a subspace of Pol<fc(5"~^) which is isomorphic to Harm^'^^. More generally, 
the set {(p{f)P{u) : f G Harm^-\degP < i} is a subspace of Pol<fc+i(5"-i) 
which is isomorphic to z + 1 copies of Harm^^^. By induction on k and i there 



8 CHRISTINE BACHOC AND FRANK VALLENTIN 

exist polynomials Pi{u) of degree i such that ip{Hj^~^)Pi{u) = Hj^^^^. Note that 
this construction could be used to derive decomposition ^ explicitly. 

We can exploit the fact that the subspaces HJ^J^ are pairwise orthogonal to prove 
an orthogonality relation between the polynomials Pj. Then this orthogonality 
relation will enable us to identify the polynomials Pj as multiples of Gegenbauer 
polynomials. Let us recall that the measures on S^~^ and on S"*^^ are related by: 

Whenever i j we have for all / S HJ^~^ 

= — / ipif)P^{u)ip{f)Pj{u)du;n{x) 

= — [ /(C)'(l - u'fn{u)P,{u)dcOn{x) 

= — I /(C)'d^n-l(C) /' (1 - V?f^^^'''^l''P,{u)P,{u)du, 

J-1 

from which we derive that 

(1 - u2)'=+("-3)/2p.(u)P, ('u)du = 0; 

hence the polynomials Pi{u) are proportional to P""^^'^(n). We obtain an orthonor- 
mal basis of Hj!~^^- from an orthonormal basis (/i, . . . , fh) of Hj^~^ by taking 

= \i^{f j)P^^'^^ {u) for a suitable normalizing factor Aj. We compute Aj in a 
similar way as above: 

1 = — / {K^{fj)Pl'^''Hu)fdu:n{x) 

^ f iMOfdu^n-iiC) C Af (1 - u'f^^-'^'\Pr'\u)fdu 



iOr. 



A2(1 _ ^x2)fc+{n-3)/2(^n+2fc(^)^2^^^ 



From the addition formula dUl applied to (P" i'^))'^ one easily shows that 



1 '^n+2fc-l"j 

SO we obtain 

l^n-l '^n+2fc * 

Now we are in the situation of Theorem [TT] with 

and their orthonormal basis (eo,i, . . . , eo,/i)) • • • > (cd-fc,!, • • • , ^d-k^h)- The isomor- 
phisms 0j are the multiplications by (Ai/Ai)P"^^'^(ii). 

Then, the coefficient with < i,j < d — k, of is given by: 
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1 ^ 

s=l 

h 

s=l 

s=l 

where we have written y = ve + Vl — v'^$, and where we apphed the addition 
formula (01) to get the last equality. Now we define Xij = XiXj and since 

C-C = {t- H/V(i-^')(i-^'), 

we have completed the proof. □ 

Remark 3.3. We would like to point out that the role of the number d is only to 
cut to a matrix of finite size. Indeed, d does not enter in the expression of 
(Yj^^ . ^{u,v,t). It is better to view the matrices Y^ as matrices of infinite size 
with all finite principal minors having the matrix-type positivity property. 

Remark 3.4. For the semidefinite programming bounds in Section E] we only use 
the matrix-type positivity property of the matrices Y^. This property is preserved if 
one replaces Y^ by AY^A^ with an invertible matrix A. So, e.g. one could replace 
the expression of (Y"^") . v, t) in (fT2l) by the simpler u^v^Q^~^{u, v, t). 

Due to the specific choice of the unit vector e defining the subgroup H, the co- 
efficients of Yj^ are not symmetric polynomials. We introduce the symmetrization 
5^ of Y^ and state the announced property 

Corollary 3.5. For all d > 0, for all k > 0, let Y^ be the matrix in Theorem \3.2\ 
and let SJ^ be defined by 

(13) 5," = i^an", 

where a runs through the group of all permutations of the variables u, v, t which 
acts on matrix coefficients in the obvious way. Then the matrices 5^ are symmetric 
and have symmetric polynomials as coefficients. We have: 

(14) For all finite C C ^ Y^{e ■c,e-c',c- a') >z 0, 
and 

(15) for all finite C C S''~\ ^ 5^(c • c', c • c", c' ■ c") t 0. 

{c,c',c")ec3 
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Proof. Note that {Y^^ .^{u, v, t) = (Yk\. (^j '^^ which gives the desired prop- 
erties of S'^. Property (ITOl ) rephrases to ([l4l ) and property ([T5] ) is obtained from 
([141 ) by taking e = d' ^ C and summing over all c" G C. □ 

To end this section we show that the positivity property ([Hi is actually a con- 
sequence of the matrix-type positivity property As shown in the following 
proposition one can express the polynomials PJ! as a linear combination of diago- 
nal elements of the matrices YJ} with non-negative coefficients. 

Proposition 3.6. We have the following expression for the polynomials PJ^it) in 
tenns of matrix coefficients ofY^{u, v, t): 

(16) n"w = E-^(^").-.,.-.(-'^'*)- 

Consequently, property (1141 ) or property ([2]) implies ([T]). 

Proof. The addition formula ((4]) holds for any orthonormal basis of HJ}. We take 
an orthonormal basis of HJ!^ obtained by concatenation of orthormal basis of the 
spaces i^ofe^ • • ■^^k'k'- (^s i,ej2> • • • denotes an orthonor- 

mal basis of -ff"^^, we have from Theorem 13. II 

(^s")fc_,,fc_,(e-x,e-y,x-y) = — ^ ^ ej,i(x)ej,i(y). 

i=i 

By the addition formula ([4| 

k h 



Pkix■y) = ^i2^2<^('^K^(y) 



un ^ ^ 

k s=0 i=l 

k > n— 1 



s=0 "'k 



and hence 



k 1 n— 1 



uri 



s=0 '^k 



Since the coefficients are non-negative, and since the diagonal ele- 

ments of a semidefinite matrix are non-negative, ([T]) is a consequence of ([T4l ). 
With the action of the permutation group of the variables u,v,t 

1 k un—l 

-^{Pk{n)+pi:{v)+pm = 

s=Q k 

Replacing u = c- d ,v = c- c" , t = c' ■ c" and summing over (c, c', c") G (7^ for a 
code C, we obtain ([B from ([2]). 

□ 
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4. The Semidehnite Programming Bound 

In this section we set up an SDP whose optimum gives an upper bound for 
A{n,9) which is at least as good as the LP method. 

For a spherical code C we define the three points distance distribution 

x{u, V, t) := card{(c, c', c") E C"^ : c ■ c' = u, c ■ c" = v , c ■ c" = t}, 

CcLrd.(^u J 

where u,v,t S [—1,1] and the matrix 




being the Gram matrix of three vectors on a unit sphere, is positive semidefinite. 

The last condition together with the first is equivalent to the fact that the deter- 
minant of the Gram matrix is non-negative, hence 

(17) 1 + 2uvt -u^ -v^ -t^ >0. 

The two point distance distribution x{u) as defined in Section[2]and the three point 
distance distribution x{u,v,t) are related by x{u,u, 1) = x{u). The three point 
distance distribution satisfies the following obvious properties: 

x(n, V, t) > 0, 
x(l,l,l) = 1, 

x{a{u), a{v), a{t)) = x{u,v,t) for all permutations o", 
x(n, V, t) = card(C)^, 

u,v,t 

x{u, u, 1) = card(C). 



Furthermore, from the positivity properties Q and ([T5] |. we have for any d >0: 

^x{u,u,l)P^{u) >0 fork = l,...,d, 

u 

^x{u,v,t)S]i{u,v,t) to foTk = 0,...,d, 



u,v,t 



where the matrix 5^' has size {d — k + I) x {d — k + 1). If the minimal angular 
distance of C is 0, we have moreover 

x{u, v,t) = whenever u,v,t ^ [—1, cos 9] U {1}. 

To factor out the action of the permutations of the variables u, v, t we introduce the 
domains 

D = {{u, v,t) : -1 < u < V < t < cos6 and 1 + 2uvt - v? - v"^ - > 0}, 
Dq = {(u, u, 1) : — 1 < u < cos 9}, I = [—1, cos 9], 
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and m{u, v, t) with 

m{u, v,t) = < 



6 if u ^ V ^ t, 

3 if u = Vy^t or Uy^v = t or u = ty^v, 
1 if u = V = t. 



From the discussion above, a solution to the following semidefinite program in the 
variables x'{u, v, t) = m{u, v, t)x{u, v, t) is an upper bound for A{n, 9): 



1+ max I - x'{u, u, 1) : 



x'{u, v,t) = for all but finitely many (u, v,t) e DU Dq, 
x'{u, v,t) >0 for all {u, v,t) e DU Dq, 

ue/ {u,v,t)&D 

3 + E ^' (^) - for /c = 1, . . . , d, 

5^(1,1,1)+ ^ x'(M,u,t)5^(n,f,t) ^0 forfc = 0,...,d}. 

(M,D,t)gDUDo 

The third constraint deserves some further explanation. We have already noticed 
that 

card(C)^ = l+ ^ t) = (^1 + ^ x(n, n, 1)^ ^, 

{u,v,t)£D\jDo u£l 

which implies 

^ x'{u,v,t) + ^Yx'{u,u,l) - (^^Yx'{u,u,l)^ > 0, 

{u,v,t)eD uei uei 

and this is equivalent to the semidefinite condition: 

/ 1 ^^x'(n,n,l) ^ 



^Yx'{u,u,l) ^ x'{u,v,t) + ^Yx'{u,u,l) 

V «e-f (u,v,t)eD u&i 



>- 0. 



Remark 4.1, We wawf ?o /7o/«f out that, despite of the fact that Q implies ([T|), a5' /i' 
proved in Proposition \3.6\ the inequalities 3 + '^,^^j x' {u,u,l)PJ^{u) > should 
not be removed from our SDR Indeed, the last inequalities do not imply them for 
an arbitrary set of numbers x'{u, v, t), unless these numbers satisfy the additional 
equalities: 

x{u, v,t) = [ x{u, u, 1)) x{t, t, 1) for all t. 



NEW UPPER BOUNDS FOR KISSING NUMBERS FROM SEMIDEFINrrE PROGRAMMING 13 



These equalities do hold for codes, but they are not semidefinite conditions. It can 
be noticed that the third constraint in the maximization problem above is a weaker 
consequence of them. 

Just like in the LP method, the main problem with the above SDP, is that the 
unknowns x{u,v, t) are indexed by a continuous domain of M^. We cannot exploit 
the information that only a finite number of them are not equal to zero, because we 
don't know to which values of (u, v, t) they correspond. We solve this problem by 
applying duality theory. 

Before we derive the SDP dual to the above one we recall the principle of weak 
duality. We use the standard notation for the inner product of symmetric matrices: 
{A, B) = Trace(v4B). Let J be a (possible infinite) set of indices, let S.j € ]R"^x™ 
be symmetric matrices with j € J, let C G ]^mxm ^ symmetric matrix, and 
let Cj G M be real numbers. Suppose that the real numbers Xj G M are a feasible 
solution of the primal problem, i.e. Xj = for all but finitely many j G J, and 
C — "^j^j XjSj >z 0. Furthermore, suppose that the symmetric matrix F G W^^"^ 
is a feasible solution of the dual problem, i.e. {F, Sj) = Cj for all j G J, and F ^ 0. 

Then, we have ZljeJ '^j^i = ^^ieJ ^i'^i' - -^)- 

In our case this specializes as follows: The set of indices is J = Dq U D. The 
matrices 5'(„ t, are block matrices with four blocks of different type. We get one 
block for each positivity constraint in the above SDP. So F is also a block matrix 
with four blocks of different type. In this case it can be simplified to three blocks. 

The first block of F consists of the matrix ( ) . The second block of F is the 

\ 012 022 J 

diagonal matrix with coefficients ai , . . . , The third block of F is again a block 
matrix with blocks Fq,. . . ,Fa which have the same size as the matrices 5^. The 
matrix C is a block matrix as well. The first block of C contains the matrix ( q o )• 
The first entry of the second block is 3, the other entries in this block are zero. The 
third block of C consists of the matrices 5*^(1, 1, 1). The real numbers C(„ j^^^) are 
zero if {u,v,t) G D and equal to 1/3 if {u,v,t) G Dq. In the following theorem 
we give the SDP dual to the above one. Furthermore we apply the simplification 
5^(1,1,1) = Ofor A; > 1. 

Theorem 4.2. Any feasible solution of the following semidefinite problem gives an 
upper bound on A{n, 9): 

d 

l + min{ ^afc + 6ii + (Fo,5?(l,l,l)) : 

k=l 

( foil 6i2 \ n 
\ bl2 b22 J 

Ofc > for k = 1, . . . , d, 
FkhO fork = 0,...,d, 

d d 

5] afcPfc'^(n) + 2612 + 622 + 3 < -1, 

k=l k=0 
d 

b22 + ^{Fk,SUu,V,t))<0}, 
k=0 
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where the last inequality holds for all {u,v,t) £ D and the second to last inequality 
holds for all u G I. 

Note that if the last inequaUty holds for all {u, v, t) G D, then it also holds for 
the larger domain 

D' := {{u, V, t) : -I <u,v,t < cos and 1 + 2uvt -v? -v"^ -t^> 0}, 
because the coefficients in are symmetric polynomials. 

5. Computational Results 

In this section, we describe one possible strategy to derive explicit upper bounds 
for r„ from Theorem 14.21 Thereby we make use of techniques form polynomial 
optimization introduced e.g. in lfT3l and ifTTl which we shall briefly recall here. 

We consider the polynomials 

p{u) = -(u + 1/4)2 + 9/16, 

Pl{u, V, t) = p{u), P2(U, V, t) = p{v), P3{U, V, t) = p{t), 

Pa{u, v,t) = 1 + 2uvt — — — t'^, 
and we obviously have 

/ = {ueR:p{u)>0}, 
D' = {{u,v,t) eR^ ■.pi{u,v,t)>0, i = l,...,4}. 

We say that a polynomial / G M[xi, . . . , Xn] is a sum of squares if it can be 
written as / = Yli=i9i' for A; e N and gi G M[xi, . . . , A polynomial 
p{xi, . . . ,Xn) of degree 2m is a sum of squares if and only if there is a posi- 
tive semidefinite matrix Q so that p{xi, . . . = z^Qz where z is the vector 
of monomials z = (1, xi, . . . , Xn, xiX2, ■ ■ ■ , Xn-iXn, . . . , x™). So assuring that a 
polynomial is a sum of squares is a semidefinite condition. 

It is easy to see that the last two conditions of the semidefinite program in The- 
orem |42] are satisfied if the following two equalities hold: 

d d 
-1 - ^ akP^iu) - 26i2 - 622 - 3 Y^i^k, S'^iu, u, 1)) = q{u) + p{u)qi{u), 

k=l k=0 

d 4 

-b22 - '^i^k, SJ^{u,v,t)) = r{u,v,t) + y^^pi{u,v,t)ri{u,v,t) 

k=0 i=l 

where q, qi and r, ri, . . . , 7-4 are sums of squares. 

It is not apriori clear that the relaxation of using this specific sum of squares 
representation is strong enough. The following theorem of M. Putinar justifies our 
approach. 

Theorem 5.1. (^) Let K = {x e W : pi{x) > 0, . . . ,ps{x) > 0} be a 
compact semialgebraic set. Suppose that there is a polynomial P of the form P = 
q + Piqi + • • • + PsQs^ where q and all qi's are sums of squares, so that the set 
{x G M" : P(x) > 0} is compact. Then, every polynomial p which is positive on 
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K can be written as p = r + piri + • • • + Psrg, where r and all Vi 's are sums of 
squares. 

Now we use these considerations to formulate a finite-dimensional semidefinite 
program which gives an upper bound on the kissing number r„: We fix d and 
restrict the polynomials q, qi,r,ri, . . . ^r^ to polynomials having degree at most 
N, with N > d. Then we can use the computer to find a feasible solution of this 
finite-dimensional semidefinite program. A feasible solution of it is at the time a 
feasible solution of the SDP in Theorem 14.21 So it gives an upper bound on the 
kissing number t„. 

We implemented this approach and give our results in Table 5.1. 





best lower 


best upper bound 


LP 


SDP 


n 


bound known 


previously known 


method 


method 


3 


12 


12 (Schutte, v.d. Waerden [23]) 


13 


12 


4 


24 


24 (Musin [15]) 


25 


24 


5 


40 


46 (Odlyzko, Sloane HM) 


46 


45 


6 


72 


82 (0.,S. L16J) 


82 


78 


7 


126 


140 (0.,S. El) 


140 


135 


8 


240 


240 (0., S. fT6l, Levenshtein fTT]) 


240 


240 


9 


306 


379 (Rzhevskii, Vsemirnov 1,21,1) 


380 


366 


10 


500 


594 (Pfender mi) 


595 


567 



Table 5.1. Bounds on r, 



The values of the last column were found by solving the above semidefinite 
program for the values d = 10 and = 10. The values of the third column 
were obtained by Odlyzko and Sloane by Theorem 12.11 using the value d = 30. 
They pointed out that even d = 11 would suffice for n < 10. Our calculations 
were performed by the program csdp (Version 5.0) of B. Borchers ([2]) which is 
available on the Internet (http://infohost.nmt.edu/~borchers/csdp.html ). After solv- 
ing the SDP with csdp we checked independently whether the solution satisfies 
the desired constraints. This can be done using rational arithmetic only. So our 
computations give rigorous proofs of the stated upper bounds. Due to numerical 
instabilities we were not able to perform this calculation for larger n and/or larger 
d, N. The smallest values of d and N which solve the kissing number problem in 
dimension 3 is d = A^ = 5. Then, we obtain by the SDP method rs < 12.8721. 
For the kissing number problem in dimension 4 it is d = A^ = 7, and the SDP 
method gives < 24.5797. 

For the lower bounds in the first column we refer to the Catalogue of Lattices of 
G. Nebe and N.J.A. Sloane (http://www.research.att.comr njas/lattices/kiss.html ). 

Using the polynomial p{u) = — (n + 1/3)^ + 4/9 we computed upper bounds 
for A{n, cos^^ 1/3). Hereby we improved several entries of the Table 9.2 of JTl 
where all best upper bounds previously known were obtained by the LP method. 
We give our results in Table 5.2. Again we used the values d = 10 and A^ = 10 to 
obtain the last column. 



16 



CHRISTINE BACHOC AND FRANK VALLENTIN 





best lower 


best upper bound 


SDP 


n 


bound known 


previously known 


method 


3 


9 


9 


9 


4 


14 


15 


15 


5 


20 


24 


23 


6 


32 


37 


35 


7 


56 


56 


56 


8 


64 


78 


74 


9 


96 


107 


99 


10 




146 


135 



Table 5.2. Bounds on A{n, cos ^ 1/3). 



We were also able to improve the best known upper bounds for the so-called 
Tammes problem with N spheres: What is the largest minimal angle 9{N) that 
can be obtained by a spherical code of with cardinality N. Let us recall that 
the answer is only known for < 12 and for = 24 (see 171 Ch. 1]). For 
N = 13, the best known lower bound is 0.997223593 57.1367031° whereas 
the best-known upper bound is 1.02746114 < 58.8691870° due to K. Boroczky 
and L. Szabo iffl. We obtained ^(3, cos-^ (0.5225)) < 12.99 using d = N = W, 
giving the new upper bound of 1.02101593 ^ 58.4999037°. Other values are col- 
lected in Table 5.3; the lower bounds are taken from the homepage of N.J.A. Sloane 
(http://www.research.att.com/~njas/packings/). The upper bounds for > 14 
where established in |5j. 





best lower 


best upper bound 


SDP 


A^ 


bound known 


previously known 


method 


13 


57.13 


58.87 


58.50 


14 


55.67 


58.00 


56.58 


15 


53.65 


55.84 


55.03 


16 


52.24 


53.92 


53.27 


17 


51.09 


52.11 


51.69 



Table 5.3. Bounds on ^(A^) (given in degrees). 
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